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Abstract 

The celebrated Leibnitz triangle has a remarkable property, namely that each of its elements 
equals the sum of its South- West and South-East neighbors. In probabilistic terms, this corre- 
sponds to a specific form of correlation of N equally probable binary variables which satisfy scale- 

invariance. Indeed, the marginal probabilities of the iV-system precisely coincide with the joint 

i_ r°° [p(x)]i 

probabilities of the (N — l)-system. On the other hand, the nonadditive entropy S q = 

(q £ K; S\ = — p(x) lnp(x)), which grounds nonextensive statistical mechanics, is, under ap- 
propriate constraints, extremized by the (q-Gaussian) distribution p q (x) oc [1 — (1 — q) (5 x 2 ] 1 /^- 1 ?) 
(q < 3; p\(x) oc e - ^ 2 ). These distributions also result, as attractors, from a generalized central 
limit theorem for random variables which have a finite generalized variance, and are correlated in 
a specific way called q-independence. In order to physically enlighten this concept, we introduce 
here three types of asymptotically scale-invariant probabilistic models with binary random vari- 
ables, namely (i) a family, characterized by an index v = 1, 2, 3, ... , unifying the Leibnitz triangle 
{u = 1) and the case of independent variables [y — > oo); (ii) two slightly different discretizations of 
g-Gaussians; (iii) a special family, characterized by the parameter x> which generalizes the usual 
case of independent variables (recovered for \ = 1/2)- Models (i) and (iii) are in fact strictly 
scale-invariant. For models (i), we analytically show that the N — > oo probability distribution is 
a g-Gaussian with q = {v — 2)/{u — 1). Models (ii) approach g-Gaussians by construction, and 
we numerically show that they do so with asymptotic scale-invar iance. Models (iii), like two other 
strictly scale- invariant models recently discussed by Hilhorst and Schehr (2007), approach instead 
limiting distributions which are not g-Gaussians. The scenario which emerges is that asymptotic 
(or even strict) scale-invariance is not sufficient but it might be necessary for having strict (or 
asymptotic) g-independence, which, in turn, mandates (/-Gaussian attractors. 

PACS numbers: 05.20.-y, 02.50.-r, 05.70.-a 

Keywords: Central limit theorem, q-Gaussians, Nonextensive statistical mechanics, Nonadditive entropy 
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I. INTRODUCTION 



The central limit theorem (CLT) provides a most powerful tool to explain the ubiquity 
of Gaussian distributions in physical systems. It explains that N independent or weakly 
correlated arbitrarily distributed random variables, with finite variances, sum up to Gaus- 
sian probability distributions for iV — > oo, corresponding to the thermodynamical limit in 
physical systems. This theorem constitutes part of the foundations of Boltzmann-Gibbs 
(BG) statistical mechanics, making it possible to describe a vast number of systems without 
accounting for the specific micro-dynamics constituting them. 

On the other hand, in systems dominated by strongly correlated microscopic events, 
correlations have remained a stumbling block for the researcher, making it extremely difficult 
to adequately take into account the contribution of all the microscopic events in order to 
get a general macroscopic behavior. 

Recently, specific generalizations of the CLT have been proposed taking into account some 
classes of global correlations, typically correlations over long distances 

flaaaaaaaa. 

Let us briefly review the present situation. When we have N identical independent random 
variables whose individual distribution has a finite variance, their sum approaches, as N 
diverges and after appropriate centering and scaling, a Gaussian distribution. This is the 
so-called standard CLT. If the individual variance diverges (due to fat tails of the power-law 
class, excepting for possible logarithmic corrections), the attractor is a Levy distribution 
(also called sometimes a-stable distribution). If the variables are not independent but q- 
independent (the q = 1 particular instance recovers standard probabilistic independence), 
;hen, if certain g-generalized variance is finite, the attractor is a q-Gaussian distribution 
2I; if that variance diverges (due to specific power-law asymptotics), then the attractor is a 
(g, a) -stable distribution [3J. These various results have been numerically illustrated in 

aa, 



and some extensions can be seen in 



101 ]. When q — 1, the correlations disappear, 



and the g-Gaussian ((g, a)-stable distribution) reproduces a Gaussian (Levy distribution). 



In terms of mathematical grounding of statistica 



nonextensive statistical mechanics 
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mechanics, the q 7^ 1 CLT cases play for 
171 ] the same role that the q = 1 CLT cases 



play for the BG theory. The development of this theory is motivated by the observation 
that g-Gaussians (or distributions very close to them) appear in many real physical systems, 
such as cold atoms in dissipative optical lattices la ], dusty plasma [lj]], motion of Hydra 



3 



cells 



20| . and defect turbulence [2l|. This suggests that this kind of probability distribution 



plays an important role in system out-of equilibrium presenting global correlations. 



A centra 
(defined in 



point of this generalized theorem is of course the hypothesis of g-independence 



2J through the g-product 



22 



23| , and the g-generalized Fourier transform ) 



This corresponds, when q ^ 1, to a global correlation of the N random variables. Its rigorous 
definition is however not transparent enough in physical terms. An important goal along 
these lines is therefore to describe in simple terms the basic physical assumptions behind 
the mathematical requirement of g-independence. Two types of simple models have been 



recently introduced JjJ, LL2| in order to provide this insight. They are hereafter referred to as 



the MTG and the TMNT models respectively. The first one is a discrete model with binary 
random variables; the second one involves continuous variables. They both are strictly scale- 
invariant and have been found numerically to converge, when N increases, on distributions 
remarkably close to g-Gaussians. However, a rigorous analytical treatment showed that 
the functional form of the probability distribution, although being amazingly well fitted by 
g-Gaussians in both models, differs from it in the thermodynamical limit 131 ] . This fact 
established that strict scale-invariance (hence asymptotic-scale invariance) is not sufficient 
for having g-Gaussians as limiting distributions. It remained open the question whether 
scale-invariance allows for such limiting distributions. In the present paper, we precisely 
clarify this central issue, thus providing some insight into this problem. 

After some brief review of the theoretical frame within which g-Gaussians emerge, we 
introduce, in Sec. a strictly scale-invariant family of Leibnitz-like triangles, for which 
the limiting distribution can be exactly obtained. Its limiting distributions are g-Gaussians, 
which proves that scale-invariance is consistent with g-Gaussianity. In Sec. IIVI we dis- 
cretize (in two slightly different manners) g-Gaussians. We then show numerically that 
these discretized distributions approach the limiting ones — g-Gaussians by construction — 
with asymptotic scale-invariance. This illustrates that both strict and asymptotic scale- 
invariances are compatible with g-Gaussianity. Finally, in Sec. |V] we introduce another 
family of strictly scale-invariant probabilistic triangles which, like the MTG and TMNT 
models, do not converge onto g-Gaussians, but onto rather curious distributions, having a 
singular behavior at infinity. We conclude in Sec. IVIl 
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II. g-GAUSSIANS 



In addition to their appearance in a g-generalized CLT, g-Gaussians are the nonextensive 



statistical mechanical 



14 
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, [l7| analog to Gaussians in the BG theory. By introducing 
a generalized entropic functional, a generalized thermostatistics could be developed that 
exhibits a thermodynamic scenario similar to that of the original one. This theory accounts 
for a class of systems where the BG theory fails. The entropy S q (with q ^ 1) was proposed 
as an alternative to the BG entropy for complex systems, e.g., trapped in nonergodic non- 



equilibrium states |26l.l27]. or nonlinear dynamical systems at the edge of chaos 28|,|29|]. Also, 
a connection between this generalized entropy and the generalized nonlinear Fokker-Planck 



30). 



equation leading to anomalous diffusion has been established 

Indeed, pi(x) oc e~ x2 ^ 2a2 optimizes the BG entropy S\ = —k J™ p±(x) \npi(x) dx, under 
constraints dxpi(x) = 1 and (x 2 )i = §°° oo dxx 2 pi(x) = a 2 . Analogously, g-Gaussians 



p q (x) oc [l 



q)(3x 2 ] 1 



(?<3) 



(1) 



optimize the entropy 



S q\Pq( X )] = k - 



dx IPgix)}* 1 



under constraints §™ dxp q {x) 



1 and [31] 



q-1 



dx X 2 [p q (x)] q 
<d 

dx \p q (x)] q 



a 



It must be noted that g-Gaussians (CQ) have compact support (\x\ ^ 1/ y/ (1 — q)(3) for q < 1 
and are defined for all x for q ^ 1. In addition, the second moment of g-Gaussians remains 
finite for q < 5/3. In the following, we will consider (3=1. The distribution P q (x) = 
[p q (x)] q / j°° dx[p q (x)] q is called escort distribution {24] and its relevance is discussed, for 



instance, in [25(. S q is nonadditive for q ^ 1 since, for two independent systems A and B, 
we easily verify (assuming k = 1) S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B). 

In complex cases, the BG entropy generally looses its extensivity, i.e., it no more (asymp- 
totically) increases linearly with the system size. In this paper, the emphasis is to consider 
g-Gaussians as being limiting probability functions characteristic for non-equilibrium states. 
We will determine the characteristic entropic index for the simple systems described in the 



following sections but do not necessarily expect them to yield extensivity of the g-entropy 
with the same value for q that they exhibit in the stationary state distribution. 

III. FIRST MODEL: A FAMILY OF LEIBNITZ-LIKE TRIANGLES 

In a probabilistic context, scale invariance will be said to (strictly) occur when, for a set 
of N random variables, the functional form of the associated marginal probabilities of the 
N- variables set coincides with the joint probabilities associated with the (N — l)-variables 
set, i.e, when 

J p N (x 1 ,x 2 , x N - U x N ) dx N = p N -i{x 1 ,x 2 , • • • , x N _i) . (2) 

This relation is always valid for independent random variables, where the joint probability 
corresponds to the product of the individual probabilities, but it is by no means necessarily 
valid for correlated ones (see, for instance, Sec. IIVI for a counter example). 

We take now the case of a set of binary independent variables, each one taking values 1 or 
with probabilities p and 1 — p respectively. For N = 2, the joint and marginal probability 
distributions are given by Table [B 



X1 \ x * 


1 







1 


2 

P 


p(l -p) 


p 





P(l -P) 


(1-p) 2 


1 — p 




P 


1 — p 


1 



TABLE I: Joint probability distribution for a set of N = 2 independent binary variables. 

Last row (column) of Table [U represents the marginal probabilities of x 2 (xi) which 
reproduce the form of the probability distribution for each single (N = 1) variable. For the 
N = 3 case, it is necessary to project a cube in the plane in order to represent the whole set 
of probabilities (Table HT|) . 

Each box of Table [III contains two probabilities. The one in brackets stands for the case 
£3 = 0, the other being for £3 = 1. Adding up the two probabilities of each box of Table HT1 
we get the corresponding box of Tabled! so scale invariance, eq. (T5J), comes up again, as it 
does when increasing N. 
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Xl \x 2 \ 


h- 1 









p 2 (l-p) 


1 








[P 2 (l-P)] 


\p(l-p) 2 ] 





P 2 (1~P) 
[P(l-P) 2 } 


p(l-p) 2 

1(1 -p) 3 ] 



TABLE II: Joint probability distribution for a set of N = 3 independent binary variables. 

It is clear that among the 2 N elementary events of the sample space, only N + 1 have 

different probabilities r N ^ n = p N ~ n (l — p) n , for n = 0, . . . , N, which, as a function of N, can 

be displayed in a triangle in the form 
(N = 0) 1 

(N = l) p 1-p 

(N = 2) p 2 p(l-p) (1-P) 2 

(N = 3) p 3 p 2 (l-p) P(l-p) 2 (1-P) 3 

The probabilities rj\r,n are the joint iV-variable probabilities. The above triangle reflects 
the aforementioned scale invariance, eq.([2]), in the sense that its coefficients satisfy the 
relation 

r N,n + r N,n+1 = r N-l,n (3) 

that is, the sum of two consecutive coefficients (marginal probabilities of a iV-system) in the 
same row yields the coefficient on top of them (joint probabilities of the (TV — l)-system). 
In other words, the corresponding marginal probabilities happen to coincide with the row 
just above, a quite remarkable property (by no means general: see Sec. II VI) . More precisely, 
we are comparing two systems: one with N elements and one with (N — 1) elements. And 
eq. ([3]) means that the probabilistic observation of the (N — l)-system coincides with the 
observation of (N — 1) particles of the A^-system. Relation ([3]) can alternatively be given as 
a rule to generate the triangle together with the starting condition rjv,o = P N for each row. 

Let us focus now on the random variable z = X\ + x 2 + • • • + Xn- It takes the val- 
ues 0, 1, . . . , N, with a degeneracy — imposed by the identical character of the A^ binary 
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subsystems — given by the binomial coefficients ( ), so the actual set of probabilities for z 

p N>n = P(z = N - n) = r^j r N>n (4) 

are to be calculated multiplying the above triangle by the Pascal triangle. It can be easily 
verified that pN )n is the binomial distribution which has, as limiting probability function 
(N — > oo), a Gaussian. 

Scale invariance condition ([3]) is the so called Leibnitz triangle rule. The Leibnitz triangle 

(N = 0) 1 



(N=l) 
(N = 2) 
(N = 3) 
(N = A) 
(N = 5) 



1 1 

i ~ 2 i 5 i 

i ~ 3 j_ « j_ ~ 3 i 

l 4 1 12 ± 12 1 4 1 

l 5 ± 20 L 30 1 20 ± 5 1 

6 30 60 60 30 6 



satisfies condition ([3]) and differs from the independent case in the definition of the starting 
condition, now given by nv,o = — - ■ Leibnitz triangle coefficients may be interpreted as 
a way to introduce correlations in the N random variables system. Probabilities are again 
calculated with (TjJ, i. e., by multiplying Leibnitz and Pascal triangles to get pN, n = 
Hence, the Leibnitz triangle rule leads to a uniform probability distribution, and so can be 
related to a g-Gaussian in the limit case q — ► — oo. 

We will now generalize the Leibnitz triangle introducing a family of scale invariant tri- 
angles rjy , v = 1,2,..., with boundary coefficients given by 



r (i) 



r (2) 
' N,0 



r (3) 



N+V 
2-3 



(iV + 2)(iV + 3) 
3-4-5 
(iV + 3)(iV + 4)(iV + 5) 



r ^,o (N + v).-.{N + 2v-V {> 
which recovers the Leibnitz triangle for v = 1. Let us emphasize that this definition leads 
to i) positive, ii) symmetric, and iii) norm preserving (in the sense that ^2 n=0 PN,n = 1> f° r 
all values of N) triangles. As a second example, the triangle for v — 2 reads 

8 



(N = 0) 
(N = l) 
(N = 2) 
(N — 3) 
(N = A) 
(N — 5) 



1 1 

JL 5 i 5 JL 

i 10 ± 5 j_ 10 i 

1 ^ _2_ 10 _3_ 10 _2_ ^ 1 

3 7 1 35 3 70 3 35 1 7 3 



28 28 140 140 28 28 



It may be shown that the coefficients of consecutive triangles of the family are related to 
each other in the way 

M _ r N+2,n+l (R , 

'N,n - • V°) 
r 2,l 

Therefore, all of them can be expressed in terms of the Leibnitz triangle and so a general 
expression for the coefficients may be obtained 

M = r N+2(v-l),n+v-l = ~ 1)! m 

' N,n m r/ ,mio/»t . n -. \ /N+2(v-V\\ " V ' / 



rgU-i [(--l)!] 2 (^ + 2,-l)(^- 1) ) 
In particular, the central elements of the triangle (n = N/2 for even N) 

» (2^-1)1 



?■ 



[(.-i)!] 2 (iv+2,-i)(^;- 1) ) 



(8) 



can be used to generate the whole triangle starting from the center instead of the side. 

We will now show that not only the Leibnitz triangle but also the rest of the triangles of 
the family yield q— Gaussians as limiting probability distribution. In fact, there is a value 
Q = ?iim( l/ )> f° r which the qn m — Gaussian corresponds to the N — ► oo probability distribution 
defined by the corresponding triangle, that is 



pz={ n yz^v qum (x). (9) 

for iV — > oo (as n = 0, 1, N, we need to define x in terms of n and N, normally through 
appropriate centering and scaling). For this purpose we will express the boundary coefficients 
rjy in an alternative way by using partial fraction decomposition 

r ( "> - {2u - 1)] yViy ( v ~ l \ 1 no) 
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On the other hand, due to the scale invariance rule, any term of the triangle can be expressed 
function of the boundary terms in the form 



' N,n — 

8=0 

Introducing Eq. (TTU]) in Eq. (11 ip yields 

v-X 



3=0 i=0 \ / 1 J 

where we have made the substitution = 

oo 

By means of the relation -^q^ = J e - ^" 1 ""'*^, together with the binomial expansion 

a o 

(a + 6) n = £ (^a^ft*, Eq. ([12]) can be cast in the form 



2=0 

OO 



u—X "p n / v 

= E a f /^£(- 1 ) B " i ( n )e- (i 

3=0 ^ i=0 

,00 , oo 

f— 1 » n / \ i/— 1 „ 

= ^af /^e-^ +v+3 >^(-l) n - < rje i * = ^af dze^ N+u+ ^ z (e z - l) n 

3=0 ^ i=0 3=0 ^ 

= J>f / (fee-^vM^-i) = / dze~^ z e~ N ^ (13) 

i=0 £ 3=0 / 

where /(z) = z — y\xx(e z — 1) with y = n/N. 

For large N, the integral in Eq. (Tl~3]) can be evaluated by using the saddle point method. 
The minimum of f(z) is located at z* = — ln(l — y), with f(z*) = —(1 — y) ln(l — y)—y\xxy 
and f"{z*) = (1 — y)/y. Therefore, in the limit N — > oo 



) 2 /"(^) 



3=0 
3=0 



i/+j+JV(l-j/)) ln(l-2/)+JV 1/ In j, 




iV(l-y) 

^^flf (l-y)^ + ^-ViH-§. ( i4) 

3=0 

Concerning the limiting probability distribution, it can be obtained from the triangle coef- 
ficients through 

7>M(y) = iVpM = JV^r« « ^(1 - y)-^H y -(^+§) r M , (15) 
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where we have made use of the Stirling approximation. Inserting now Eq. (fl4|) in Eq. ([15]) 
yields 



u-l 



J>i (l-y) 1 ^'- 1 (16) 

J'=0 

The largest exponent of y in the distribution ffTol) is 2(V — 1). Hence, comparing with 
Eq. ([1]), the value of gn m for the g-Gaussian limiting distribution function can be obtained 
doing 1/ (1 - g Um ) = v - 1, i.e. 

z/-2 

gum = p (17) 



which implies a width of the compact support given by A = 2a/1/(1 — q\ im ) = 1\/v — 1. 
Equation ffTT|) can be re-written as q\ ira = 1 — — 1), which reminds many analogous 
relations existing in the literature, such as (/entropy = 1 — 1/d in l[ Q, 33]. 

The variable y is defined between and 1. However, we are interested in a centered 
distribution function defined within [—A/2, A/2] and thus apply the transformation x = 
2V v — l(y — 1/2). The limiting function now results 

1 1 " -1 a (v) f <r \ v +i~ 1 



= (2 "~ 1)! ri - ^ fl8l 

2 2 - i v^T[(^-i)!] 2 \ iz-iy ' 1 ; 

which exactly coincides with a g— Gaussian with q = gn m . 

In addition, V^(x) transforms into the Gaussian distribution for v — > oo, so we recover 
the statistical independence case. In fact, it can be verified that taking limit in Eq. (J7J) one 
gets lim^oo rjy n = 2"^, hence lim^oo = (^)2 _Ar , so the corresponding triangle is the 
one given at the beginning of this section with p — 1/2. 

Figure [T](a) shows V <yU \x) = ^4= py n as compared to the corresponding g lim — Gaussians 
for N = 500. It is apparent that the approximation becomes poorer when increasing v. 
Figure [0(b) shows the validity of Eq. (jHJ) for v = 5, the corresponding qu m = 3/4 and 
N =100, 200, 500 and 1000. Curves overlap with the q— Gaussian for greater values of N. 
The convergence is thus evident. 

Figure [2] shows that, in what concerns extensivity, the family of triangles (0) follows the 
Boltzmann-Gibbs prescription, that is, the value of q that makes the q— entropy extensive 

11 




-3-2-10123 -3-2-10123 

X X 



FIG. 1: a) Dotted lines: Probability distributions V^ u \x) = 2 J^ ^ p^ n , with p^ n given in ([9|) 
for N = 500 and v = 2, 3, 5 and 10. Solid lines: Corresponding (/ii m -Gaussians with qn m given in 
(fT7|) for the respective values qy m = 0, 1/2,3/4 and 8/9. The limiting case v — > oo (Gaussian) is 
also depicted for comparison, b) P (5) (x) for N = 100,200,500 and 1000 and the corresponding 
?lim— Gaussian with qn m = 3/4. 

is g e nt = 1 for all values of v. 

As a last remark, let us associate with the N random variables the variables <X; = 2xi — 1 = 
±1 (i = 1, 2, N), so that (cr^) = 0, VI We can straightforwardly prove that (<JiCTj) = ^pj, 
Vz 7^ j, ViV. In the limit v — > oo we recover (<Ji<Jj) = 0, as expected for independent 
variables, where no correlation exists. 




N N N 

FIG. 2: The q— entropy S q = — q " l - — as a function of iV for v = 1 (left), = 2 (center), 

and z/ = 5 (right), and typical values of g. In all cases q en t = 1. 
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IV. SECOND MODEL: DISCRETIZED g-GAUSSIANS 



We will now introduce another probabilistic model in which we impose a priori the con- 
dition that the N — > oo limit for the probability distribution is a g— Gaussian, with the aim 
to study whether (strict or asymptotic) scale invariance is also obtained. In order to verify 
if the concept of g-independence, i.e. correlations leading to g-Gaussians, can be related to 
scale invariance in probabilistic terms, relation ([3]) is expected to be satisfied at least in the 
limit N — » oo, i.e., asymptotically. 

Considering again the set of TV equally probable binary variables, the correlations will 
now be given in the form 

r N>n - — (19) 

' n=0 

where N + 1 equally spaced points in the support of the g-Gaussian p q (x), to be 

specified later. For the set of probabilities we again write pjv> = (^)rjv jn , which provides 
us with a discrete probability distribution which, by construction, follows the shape of the 
g-Gaussian p q (x). 

Concerning the way to choose the points xn,u in {ED a distinction must be made between 
cases q < 1 and q ^ 1. 

As mentioned before, for q < 1, g-Gaussians have compact and symmetric support of 
width A = 2/V(l - q). For this case, we will consider two different ways to choose the 
points x 7v,n in the support of the g-Gaussian: 

1) N + 2 discretization (Dl): In this implementation, we take xjy,n = ^min + (n + l)h, 
for n = 0, 1, . . . , N, with x m i n = —A/2 and h = A/(N + 2), i.e., explicitly, xn,u = a^min + 
h, x min + 2h, x min + 3h, x min + (N+l)h. 

2) N + 1 discretization (D2): Now, the points chosen differently: the same 
initial interval A breaks now into N + 1 equal subintervals (not iV + 2 as before) of width 
h = A/(iV + 1) and we take the values of p q (xN, n ) in the center of each subinterval, i.e. 
%N,n = %min + (2n + l)h/2. The whole set reads xn,u = Xmin + h/2,x m i n + 3h/2,x m i n + 
5h/2,... 1 x min + (2N +l)h/2. 

In contrast, for g ^ 1, the support for g-Gaussians is the whole real axis and we must 
take this into account in the fit. We will take an increasing width for the fit interval in the 
form Ajy = S(N + l) 7 , 5 being some initial width, and 7, with ^ 7 ^ 1, a parameter 
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determining the growth of the interval width (for 7 = we recover the q < 1 case). Now 
xn,u = £minjv + {n + l)h N , with 2 min ,7v = ~A N /2 and h N = A N / (N + 2). 

Despite the fact that different discretizations yield different triangles (1191) for a given 
value of q, let us emphasize that the corresponding limiting distributions PN,n tend to the 
same q— Gaussian p q (x) in the limit N — > 00. 

Now the following question arises. Do the triangles (Tl9|) satisfy relation ([3]) as the triangles 
([5]) from Sec. II III do? In other words, can g-Gaussians be related to strictly scale-invariant 
distributions? Strictly speaking, they are not, since relation ([3]) is not exactly fulfilled (except 
for the case q = with the first discretization Dl, as we will show later), but we will show 
(analytically in some cases, numerically in others) that these triangles are asymptotically 
scale- invariant, that is, relation ([3]) is satisfied for iV — > 00, or, alternatively, the ratio 

Qn,u = (20) 

fN+l,n + r N+l,n+l 

tends to 1 (or equivalently Qn,u — 1 tends to 0) as N increases. Note that Qo,o = Qi,o = 
Qi t i = 1 for arbitrary values of r ,o, T\ )fl and r 2j „. 

A. The q < 1 case 

Figure [3] shows Qn,u — 1 as a function of n for iV = 500 and different values of q, for both 
the Dl and D2-discretizations. It is clearly observed the proximity of Qn,u to 1, which is 
more noticeable in the center of the triangle. 

Quite remarkably, for q = 0, strict scale invariance is obtained in the first discretization, 
that is, Qn,u — 1 for all iV and n. This is so because, in this case, it can be proved that 
triangle ( fl9|) exactly coincides with the Leibnitz-like triangle of the family (JHJ) with v = 2, 
with associated probabilities p%] n = Qr$ n = 6[1 + n{N - n)/(N + 1)]/(N + 2)(JV + 3). 

An exact expression for r^ >n and hence Qn,u can also be obtained for the D2 discretization 
and q = 0, the probabilities being in this casepjv.n = 3[(2n + l)(2iV-2n + l)]/[(iV + l)(2iV 2 + 
AN + 3)]. Of particular interest are the central value, Q c = Qn,n/2, and the boundary one, 
Qo = Qn,o, of quotient (l20i) . being given by 

(N + 2) 2 iV(2iV 2 + 8N + 9) , 

Qc = 7 ^ — — t-- — r for odd N, (21) 

Vc (2N 2 + AN + 3)(N 3 + 6N 2 + 10N + 2) ' V ; 

(2iV + l)(iV + 2)(2iV 2 + 8iV + 9) 

^° (2iV 2 + 4iV + 3)(2iV 2 + lliV + 6) ' 1 J 
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s 



-0.0001 




FIG. 3: Q N , n - 1 

as a function of x n n for N = 500 and different values of q = -1, -1/2, 0, 1/4, 1/2 
for discretizations Dl (top) and D2 (bottom). Strict scale invariance is observed for q = and 
discretization Dl. In the rest of the cases deviation from zero is small and, as can be seen by 
results not shown here, decreases when increasing N. 
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From Eqs. ([21]) and Q22D results Q c -1 ~ A^~ 2 and Q - 1 TV -1 , respectively. Though 

this equations are only valid for q = 0, this trend is observed for any value of q < 1. 
Figure H] shows in a log-log plot Q c — 1 as a function of N for different values of q < 1 and 
both discretizations. It is clear that the decay follows a 1/N 2 power-law for large A" and 
any value of q. No substantial differences are observed between both discretizations. 

Analogously, Fig. [5] shows Qo — 1 as a function of A" and different values of q. We observe 
now a 1/N power-law. We found that this 1/N power-law transforms into the 1/N 2 when 
we do not take into account the complete interval of the compact support of the g-Gaussian 
under consideration (not shown). 

Concerning the extensivity of S q , the same behavior as in the previous systems is found. 
We get g ent = 1 no matter the value of q ^ 1 of the discretized g-Gaussian (let us remind that 
there is no reason for the value of q of the discretized g— Gaussian be equal, or even simply 
related, to the index g ent corresponding to the extensivity of the entropy S q (N)), and no 
matter the type of discretization (Dl or D2). Figure [6] shows the q— entropy for discretized 
g-Gaussians for typical values of q. The results are independent of the discretization. 

B. The q ^ 1 case 

A similar trend is observed for q ^ 1. Quotients Qo and Q c tend to 1 as A" increases 
for all values of 7, which, as mentioned before, determines the growth of the interval where 
the q— Gaussian is evaluated. In the case of the Gaussian, i.e. q — 1, it is known that 
7 = 1/2. Figure [7] shows the decay of the central quotient for q = 3/2 and different values 
of 7. Apparently, 7 = 1/2 provides the appropriate growth of the interval for g-Gaussians 
with g > 1 as well. For 7 < 1/2, one observes the power-law behavior only over some range 
whereas for 7 > 1/2 the decay follows a power-law with an exponent larger than —1. The 
boundary ratio displays the same dependence on 7. 

For 7 = 1/2 we verify a 1/N power-law. Figure [8] shows Q c — 1 and Qo — 1 for typical 
values of q and 7 = 1/2. 

In what concerns the extensivity of the entropy S q , the value of g ent remains 1 and is 
independent of 7. Figure [9] shows the S q (N) for typical values of q ^ 1. 



1(3 





FIG. 4: Central ratio Q c — 1 as a function of N for discretized q— Gaussians with q = — 1, 0, 1/4, 1/2 
and 3/4 for discretizations Dl (top) and D2 (bottom). Exact result (|2ip for g = is also shown 
for D2 discretization. The power law with exponent —2 is shown for comparison. 
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N 



FIG. 5: Boundary ratio Qq — 1 as a function of N for discretized q— Gaussians with q = 
-1,0,1/4,1/2 and 3/4 for discretizations Dl (top) and D2 (bottom). Exact result ([22]) for q = 
is also shown for D2 discretization. The power law with exponent —1 is shown for comparison. 
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N N N 



FIG. 6: q— entropy S q = " = q _" — , with r^ >n given in (fl~9l) . as a function of N for discretized 

g-Gaussians with q = —1 (left), q = —1/2 (center), and q = 1/2 (right) and discretization Dl. 
Results with discretization D2 are indistinguishable. In all cases q en t = 1. 




N 



FIG. 7: Central ratio \Q C — 1| as a function of N for g = 3/2, 5 = 2 and 7 =0.4, 0.5 and 0.6. 
Asymptotic power law behavior with exponent dependent on 7 is observed for 7 ^ 1/2. For 7 = 1/2 
Q c — 1 ~ iV -1 , the decay being slower for greater 7. 

V. THIRD MODEL: ANOTHER FAMILY OF GENERALIZED TRIANGLES 

As seen in Sections lIIII and HVt strictly as well as asymptotically scale-invariant probab ility 



models may lead to q— Gaussian limiting distributions. Nevertheless, as we already now 



19 





FIG. 8: Central quotient (top) and boundary quotient (bottom) decay as a function of N for 
discretized q— Gaussians with q = 1.01, 1.25, 1.25 and 2, 7 = 1/2 and 5 = 2. In all cases the trend 
Qo — 1 ~ Q c — 1 ~ iV" 1 is observed. 

scale-invariance does not guarantee q— Gaussianity. In this Section, we present a last model 
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9=0.6 



9=0.8/ 




9=0.6 



9=0.8/ 





FIG. 9: q— entropy as in FIG. E]for discretized q— Gaussians with q = 1 (left), q = 3/2 (center), 
and q = 2 (right). In all cases q ent = 1. 



to emphasize this point. 

The following strictly scale-invariant triangle 
(N — 0) 1 

1 1 

(N = l) 
(N = 2) 
(N = 3) 
(N = A) 
(N = 5) 



± 2 ± 2 £ 

£ 14 £ 14 £ 14 A 

£ 28 £ 28 ± 28 £ 28 L 

11 56 3 56 3 56 3 56 3 56 11 
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112 



112 



112 



112 



with coefficients (x = 3/7) given by 



il-iV 



X2 



1-iV. 



); n = 0,N 



(23) 



corresponds to a different way to introduce correlations in the system. In order to get 
nonnegative coefficients, the parameter x is kept within [0, 1/2]. 
The probabilities are given by 



A/ 



X2 



l-JV 



(24) 



The case of binary random variables [y — > oo of the triangles analyzed in Sec. IHIj) is 
reproduced here for x = 1/2, hence pN, n = ( n )2 _JV • 

In order to calculate the limiting probability function, the CLT states that the new 
variable x = (n — N/2)/(\/~N /2) provides a normal distribution in the limit N —>■ oo for 
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the second term of Eq. (T24|) . In addition, two delta peaks appear after substitutions 5 n fi — > 
6(&x + % ) = ^5(x + VN) and -> 5(^-f ) = ->=5{x-^N). Finally, by taking 

limits in Eq. (1241) . we obtain 



P(a;) = lim ^ P 7v,n = 2x^e"" 2 / 2 + (\ - x] lim - v 7 ^) + £(z + v 7 ^)) (25) 

JV— >oo 2 V27T / 

which consists of a Gaussian distribution plus the additional contribution of the delta peaks 
corresponding to a concentration of probability on the two sides of the triangle. 

As in the previous sections, the BG-entropy is extensive for this triangle as well. This 
may be proved directly by inserting coefficients (I24p into 

Si = - e ( ) rN > n ln rN > n ( 26 ) 

n=0 ^ ' 

yielding 

N-l 



1 /N\ 
Si = -(1 - X) ln[(- - X ) + X2'- N ] ~ X^~ N E ( n ) Ntf 1 " 

n=l ^ ' 



AT 



)] 



= (X - 1) M(i - X) + tf 1 "") - 2 X (ln(2 X ) - 2^ In^ 1 "")) + 2 X iVln2 

oc iV (27) 

for large N. 



VI. CONCLUSIONS 



A family of Leibnitz-like triangles, leading to g-Gaussians as limiting probability dis- 
tribution functions with q ^ 1, was introduced, where the limiting distribution could be 
exactly calculated. These systems correspond to iV correlated binary random variables, the 
index q characterizing the strength of correlation. The case q — > — oo corresponds to very 
strongly correlated variables giving a uniform limiting distribution. 

On the other hand, the coefficients of another type of triangles were constructed by 
discretizing g-Gaussians. These triangles, having now by construction g-Gaussians with 
q < 1 as limiting probability functions, showed a behavior with depends on the specific 
discretization of the support interval. Except for one particular case, the Leibnitz rule, 
related to system size scale invariance of the probabilities, is only asymptotically satisfied. 
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The system approaches scale invariance with a l/N 2 power- law for large N, except for 
the boundary coefficients where the convergence to scale invariance is much slower, of the 
type l/N. The l/N 2 law makes a crossover into a l/N one over the entire triangle when 
considering g-Gaussians with q ^ 1. 

Finally, another family of strictly scale-invariant triangles with a rather strange limiting 
distribution function was introduced. In the limit N — > oo, the triangles yield a Gaussian 
distribution together with two delta peaks centered at points going to infinity. 

The BG-entropy remains extensive for all three types of triangles, equally to previously 
studied Leibnitz-like triangles This may be the result of the simplicity of the models 
presented in this paper. More sophisticated models, as for instance the Hamiltonian mean 



field model (see for instance ref. 
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271]). appear to approach a g-Gaussian characterized 



by a non-equilibrium stationary state with the g-entropy possibly being extensive for q ^ 1. 
However, in the present effort we are here not particularly interested in the general relation 
between the extensivity of the entropy and stationary-state probability distributions, but 
we rather searched to find out which kind of correlation between the microscopic events of 
a system leads to g-Gaussians as limiting distributions (possibly, as attractors). 

The Leibnitz rule provides a simple tool to study models composed of correlated binary 



random variab 
addressed in 



es, and enabled the exact calculation of their limiting functions. As already 



131 ]. this rule cannot be uniquely related to nonextensive thermostatistics. 



Indeed, Leibnitz-like triangles exist which precisely lead to g-Gaussians (as shown in the 
present paper) as well as to other limiting probability functions (as shown in [l^], and also 
here) . Additionally, the present second family of triangles (with asymptotic but not strict 
scale invariance) also tended to a g-Gaussian. The scenario which emerges is that asymptotic 
validity of the Leibnitz rule might represent a necessary but surely not sufficient condition 
for the system to tend to g-Gaussians as limiting distributions when N — > oo. 

The fact that different implementations of correlations between the variables of a system 
can lead to the same function, — g-Gaussians in the present case — , can be seen as a hint 
for these functions being attractors for a variety of different systems, and so supports the 

nnnn 

demand of generality of the g-generalized central limit theorem presented in |2J, y, |4|. 
However, to assure the applicability of this central limit theorem, the stability of the q- 
Gaussians as limiting functions of the systems presented here needs to be proved, either 
by establishing that the correlations correspond to g-independence, or by introducing, for 
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example, weak perturbations. 
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